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Unclassified 


FOREWORD 


This  is  the  third  report  in  the  Data  Reduction  Program  Documenta¬ 
tion  series.  It  is  dated  according  to  the  date  of  completion  of  the  documentation.  No 
implication  is  made  that  this  program  will  not  subsequently  be  modified,  amended,  or 
superseded:  on  the  contrary,  the  history  of  radar  data  processing  is  one  of  continuous 
evolution  of  techniques,  and  it  is  unrealistic  to  assume  that  steady -state  has  been 
reached.  The  PA -229  series  is  being  published  for  the  convenience  of  interested 
parties,  and  Lincoln  assumes  no  responsibility  for  the  correctness  of  the  information 
presented,  nor  for  its  currency. 

The  preparation  of  reports  in  this  series  is  under  the  Editorship  of  Charles  R. 
Bemdtson  of  Lincoln,  and  of  D.  Nessman  and  R.  French  of  Philco-Ford  Corporation. 
Inquiries,  suggestions,  corrections,  criticisms,  and  requests  for  additional  copies 
should  be  directed  to  C.  R.  Bemdtson. 


The  principal  contributor  to  this  report  was  ].  R.  Cornelius  (Philco-Ford). 

Due  to  the  intricate,  evolutionary  manner  in  which  the  programs  came  into  being,  the 
editors  regret  that  it  is  in  general  impossible  to  give  due  credit  to  all  --  mathematicians 
or  radar  analysts  or  programmers  --  who  contributed  to  the  definition  and  writing  of 
the  programs. 
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ALREAD 


I.  PURPOSE  AND  UTILIZATION 

A.  Source  of  Data 
ALTA1R1 


B.  Data  Input 

ALTAlR  Transcription  Tape 
0.  ^.^Description 

X,  ...  ^ 


The  ALTAIR  transcription  tape  is  written  in  a  variable  record  length 
format.  Subroutine  ALREAD  is  designed  to  process  transcription  tape  records  and 
extract  data  necessary  to  run  several  programs.  It  checks  for  errors,  performs 
computations,  and  stores  the  values  for  the  calling  program. 


ALREAD  is  called  by  the  following  program/fKALTOAK  (PA-229-2); 
ALTAP  (PA-229-4);  ALT10  (PA-229- 5);  ALTFENCE  (PA-229-6)  ;mhi  ALTMOP 
(PA-229-  ).  rhis  report  applies  to  these  five  programs,  except  for  the  program 
listing  given  in  Appendix  A.  This  appendix  applies  to  ALTAP,  ALT10,  and  ALTFENCE 
The  reports  for  ALTOAK  and  ALTMOP  contain  listings  for  the  version  of  ALREAD  us¬ 
ed  by  each. 


D.  Output 

4.  A  listing  of  error  messages. 

2.  A  listing  of  format  and  calibration  records  found  by  Subroutine 
THEAD  and  stored  in  common. 

3.  Additional  output.  See  calling  programs. 
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II. 


DESCRIPTION 


ALREAD  processes  data  from  the  ALTAIR  transcription  tape,  one  physical 

record  at  a  time.  The  data  is  retrieved  from  the  tape  by  a  call  to  BREAD  an  assem  - 

2 

bier  language  routine  which  reads  the  tape  at  a  rate  of  one  record  per  call.  ALREAD 
also  calls  THEAD,  a  Fortran  routine  which,  through  BREAD,  reads  the  calibration 
records  from  the  tape,  initializes  data,  and  sets  up  required  format  tables.  ALREAD 
will  return  T,  Alt,  R,*  CADJ,  Az**  and  El**  for  a  maximum  of  300  consecutive 
pulses.  It  also  returns  RCS  data  for  each  pulse  and  gate. 

A  call  is  made  to  THEAD  and  data  initialization  takes  place.  Vital  format 
tables  are  tested  and,  if  any  are  missing,  an  error  message  is  printed  and  a  return  to 
the  calling  program  made.  If  all  tables  are  found,  the  lift-off  time^  is  "npacked  and 
stored.  Successive  calls  to  BREAD  are  made  until  the  record  containing  the  start  time 
is  found.  This  record  is  unpacked  and  vital  data  are  tested  and  stored  for  use  in  com¬ 
putations  of  T,  CADJ,  RCS,  or  for  direct  passage  to  the  calling  program. 

Data  gathering  begins  at  the  first  pulse  whose  time  ^  the  start  time.  It  stops 
at  the  first  pulse  whose  time  sthe  stop  time,  or  when  data  for  the  maximum  number  of 
pulses  has  been  gathered,  whichever  occurs  first.  If  a  maximum  of  points  are  collected 
before  TPLS  s  the  stop  time,  a  code  indicates  to  the  calling  program  that  the  stop  time 
conditiou  is  unsatisfied.  The  calling  program  will  process  this  data  before  recalling 
ALREAD. 
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Pulse  times  (TPLS)  are  computed  by  the  formula: 


TPLS  =Trec  +  (Pulse  n  -1)/PRF 
where 

^REC  *S  t^ie  c*rne  c*le  ^rst  pulse  in  the  major  cycle 
Pulse  n  is  the  nth  pulse  in  the  major  cycle. 

RCS  is  computed  as  follows: 

RCS  =  Amplitude?"  +  CADJ 
where 

CADJ  =  Calibration  constant"*  -  Power  -  10  *  (sensitivity  bit) +  + 
attenuation  +  40  *  (log  R) 

Note  that  Power  and  R  are  used  in  counts  in  this  computation. 


‘"Amplitude  is  obtained  by  indexing  the  appropriate  calibration  record  (FMAMPl 
FMAMP6)!  with  the  amplitude  values  obtained  in  FMTRSP,  Item  1. 

**Found  in  Calibration  Record  FMXSEC.  ^ 

Found  in  FMTRAMA,  Item  6. 1 


to 
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III. 


OPERATION 


A.  Input 

Call  statements  for  ALREAD,  since  they  vary  from  program  to  pro¬ 
gram,  are  given  in  the  programs  which  utilize  it. 

Definitions  of  input  parameters: 


TSTART 

Start  time  of  processing  (GMT  total  seconds) 

TSTOP 

End  time  of  processing  (GMT  total  seconds) 

INTARG 

Target  no.  to  be  processed 

INPAT* 

Sampling  pattern  in  which  initial  gate  is  located 

IPOL 

Data  channel:  1  =  LC;  2  =  RC; 

3  =  Az  error;  4  =  El  error 

NOPHA 

Type  of  data:  1  =  RCS  only;  2  =  Phase  and  RCS 

NRG 

Number  of  range  gates  to  be  processed 

ISTGAT** 

Location  within  INPAT  of  initial  gate 

Input  and  Output  Parameters 

NPTS 

Output:  Number  of  pulses  of  data  returned 

Input:  Must  be  initialized  by  calling  program 
before  each  call  to  ALREAD  ^ 

1 1 

NEWPAS 

Cycle  and  error  pointer  (see  error  returns  and 
switch  settings) 

Output 

1 .  Output  Parameters 

TLIFT 

Lift-off  time  (GMT  total  seconds) 

"  Also  called  IPAT. 

**Also  called  ISG  or  ING. 

^Set  to  zero  for  first  call.  Set  to  number  of  saved  points  for  all  subsec uent  calls. 
+  +  Also  called  IAGAIN. 
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DFPG  Frequency  and  Polarization  (e.g.  VHF  LC) 

NFP  Frequency  code:  1  =  VHF;  2  =  UHF 

Stored  in  Common 


TIMES 

XSPHA 

RANGKM 

ALSAV* 

ELSAV 

IRGA 

NFPG 

CADJ*  * 

AZI 

ELE 

FSIG 

IPHASE 

DOPM 

PRD 

IPOL 

NPOL 


Pulse  times  (GMT  total  seconds) 

RCS  and  Phase  for  each  pulse  and  gate 
R 

Alt 

El  (rad) 

Range  gate  array  associated  with  XSPHA 
Frequency  code:  1  =  VHF;  2  =  UHF 
Adjusted  calibration  constant 
Az  (rad) 

El  (rad) 

£ 

RCS  (m)  for  each  pulse  and  gate 
Phase  for  each  pulse  and  gate 
Range  rate  (km/s) 

PRF 

Data  channels  wanted 
Number  of  data  channels  used 


2.  Error  Returns  and/or  Switch  Settings 


NEWPAS  =  0  Indicates  that  re-entry  to  ALREAD  is  not  desirable 
either  because  TSTOP  has  been  reached  or  because 
a  fatal  error  has  occurred.  (See  below) 


a.  "THEAD  has  defaulted  -  run  aborted" 

b.  "At  time  =  (12345.  o789)  Target  (00)  is  not  on  the  tape, 
targets  available  are  (00,  00,  00,  00,  00)" 

c.  "Either  Polarization  or  pattern  chosen  is  not  available" 


♦Valid  only  for  first  pulse  of  minor  cycle.  it  is  repeated  for  subsequent  pulses. 
**Also  called  CALOUT. 
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NEWPAS  =  1 


Indicates  a  tape  has  been  rewound. 


NEWPAS  =  44  Indicates  that  an  ALTAIR  data  end -of -file 

record  has  been  read  from  the  tape.  A  message 
is  printed  by  ALREAD.  "End -of -file  encoun¬ 
tered  on  input  tape  after  time  <12345.  6789)" 

NEWPAS  =  55  Indicates  that  a  vital  format  table  is  missing; 

processing  must  stop.  A  message  identifying 
the  table  is  printed.  "Format  table  (NAME)  was 
not  found  -  run  aborted" 


NEWPAS  =  66 


Indicates  that  an  amplitude  or  phase  table  for  a 
requested  polarization  is  missing.  A  message 
identifying  the  table  is  printed. 


NEWPAS  =  99 


Indicates  that  the  maximum  number  of  points  per 
call  have  been  stored  but  that  TSTOP  has  not 
been  reached. 


A LREAD -generated  error  messages  and  accompanying  switch 
setting. 


t 
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At  time  (12345.6789),  minor  cycle  (0),  something  is 
screwed  up  NBP*  =  (01234).  It  should  be  (43210)" 
NEWPAS  is  not  set.  ALREAD  will  continue  by  skipping 
data  in  current  record,  and  calling  BREAD  for  next 
record. 

"NPLS  =  (00)  Minor  cycle  (0),  (0)  minor  cycles" 

A  test  for  a  value  of  NPLS  ^0.  Same  action  as  in  a. 

"At  time  (12345.  6789)  R  =  (12345.678900)" 

A  test  for  a  value  of  R  ^0.  If  this  occurs,  R  is  set  to 
1.0  and  processing  continues. 


NBP  is  number  of  bits  per  pulse  and  is  coded  on  the  tape.  It  is  also  computed. 

If  the  two  values  are  not  identical,  some  pointers  are  wrong  and  the  data  may  be 
faulty.  J 
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IV. 


PROGRAM  LIMITATIONS 


NRG 

NPTS  per  call 


120  (30  for  ALTMOP  and  ALTOAK) 

^  300  for  RCS  data  (NOPHA  =  1) 
s  150  for  RCS  and  phase  data  (NOPHA  =  2) 
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V. 


PROGRAMMING 


A.  ALREAD 

A  program  listing  and  a  flow  diagram  for  ALREAD  are  given  in  Ap¬ 
pendices  A  and  B. 

ALREAD  calls  the  following  subroutines: 

B.  THEAD  (see  Appendices  C  and  D). 

THEAD  is  used  to  process  calibration  and  format  records  which 
are  recorded  before  the  data  records  on  the  ALTAIR  transcription  tape.  THEAD  lists 
format  and  calibration  records  named  in  the  common  statement.  These  are  stored  and 
unpacked  for  later  use  by  the  main  processing  program.  A  sample  THEAD  output  is 
given  in  Appendix  E. 


The  minimum  size  needed  for  the  item  array  may  be  calculated 
by  the  following^equation: 

Item  size  =  6  *  (Total  number  of  items  stored)  +  (6  *  130) 


Calling  Sequence: 


Call  THEAD  (NEWPAS,  *) 

NEWPAS  =  Tape  processing  status;  integer  variable. 

0  =  Initial  pass  thru  tape;  key  calibration  and  format 

records  stored  and  unpacked. 

1  =  Tape  has  been  rewound;  skip  calibration  and 

format  records  but  call  BREADS  to  reinitialize 
parameters  if  data  record  processing  is  to 
continue  in  main  program. 


*  =  A  return  point  specified  by  a  statement  number 

in  the  calling  program.  Used  for  aborting  job 
by  main  program  if  wrong  tape  is  mounted. 

THEAD  calls  the  following  subroutines: 


BREADS  (entries  BREADS  and  BREAD);  HDRR  (entries  HDRR  and 
NAMET);  and  FORM.^  THEAD  also  calls  WHICHV,  a  360-system  subroutine 
indicating  whether  a  job  is  being  run  under  the  time -sharing  (CMS)  or  Batch  (OS) 
systems. 
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C.  GMTUPK  (sc^  Appendix  F) 

GMTUPK  is  a  routine  to  unpack  GMT  from  the  tape.  The  call  is 

GMTUPK  (TIME). 

GMTUPK  input,  which  is  stored  in  common,  is  two  packed  data 
words,  one  containing  h,  min,  s,  and  the  other  word  containing  decimal  fractions  of  a 
second  (GMT). 

The  output  (TIME)  is  GMT  total  seconds,  ,xnd  is  contained  in  the  call 

statement. 

D.  GET  and  IGET 

2 

GET  and  IGET  are  entries  to  subroutine  GETS. 

These  routines  will  locate  any  data  item,  unpack  it,  and  interpret 
it  according  to  the  information  in  the  format  table.  They  will  return  the  item  as  a 
binary  integer  (in  the  case  of  IGET)  or  as  a  floating  point  number  (in  the  case  of  GET). 

GET  (or  IGET)  requires  three  arguments: 

GET  (Format,  Base,  Item) 

Format  Relevant  format  table  address 

Base  Base  address  of  data  block  desired 

Item  Specific  item  number 

E.  BREAD 

An  entry  to  subroutine  BREADS  used  to  read  data  records  from 

2 

the  ALT  AIR  tapes. 


REFERENCES 

1.  "ALTAIR  Data  User*  s  Manual",  LM-97,  Lincoln  Laboratory,  M.I.T. 
(to  be  published),  UNCLASSIFIED. 

2.  "Data  Reduction  Program  Documentation,  ALTAIR  Tape  Read  Package, 
(Effective:  April  1970)",  PA-229-1,  Lincoln  Laboratory,  M.I.T. 

(17  Mar^h  1971),  UNCLASSIFIED. 
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COMMON  SYMBOLS  AND  ABBREVIATIONS 


(The  units  given  for  certain  quantities  are  the  units  commonly  used  for  those 
quantities,  unless  otherwise  noted. ) 


ADT 

Alt 

APS 

ARS 

Avg 

Az 


ALCOR  Data  Tape 
Altitude  (km) 

Average  Pulse  Shape 
ALTAIR  Recording  System 
Average,  Averaging 
Azimuth  (deg) 


CADJ 

C-band 


Adjusted  Calibration  Constant  (db) 

ALCOR  frequency,  5664  MHz  (NB)  and  5667  MHz  (WB) 


El 

EOF 


Elevation  (deg) 
End  of  File 


GMT 


Greenwich  Mean  Time 


h 

Hz 


Hours 

Hertz 


in 


Inches 


LC 


Left  Circular  Polarization 


min 


Minutes 


NB  Narrow  Band 

NRTPOD  Non-real  Time  Precision  Orbit  Determination  Program 


POD 

Phase 

PRF 

PRI 

pps 

pts 


Project  PRESS  Operation  and  Data  Summary  Report 
Presented  in  deg 

Pulse  Repetition  Frequency  (pps) 

Pulse  Repetition  Interval  (s) 

Pulses  per  second 
Points 
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R 

R 

rad 

RC 

RCS 


Range  (km) 

Range  Rate  (km/s) 

Radians 

Right  Circular  Polarization 
Radar  Cross  Section  (dbsm) 


s 

SD 

w 

T 

TAL 

UHF 


Seconds 

Standard  Deviation  of  Wake  Velocity 
Time 

Time  After  Launch  (s) 

ALTAIR  Frequency;  415  MHz 


V 

vd 

V 

w 

VHF 


Velocity 

Doppler  Velocity 

Mean  Wake  Velocity 

ALTAIR  Frequency;  155.  5  MHz 


WB 


Wide  Band 


/V 

* 


Total  Off-axis  Angle  (deg) 
Wavelength 

Denotes  Multiplication 
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FLOW  DIAGRAM  SYMPOI.S 


PROCFSS.  ANNOTATION 


CD 


TERMINATOR 


SUBROUTINE:  where  NAME  is  the  entry 
call  Into  the  subroutine 


CONNECTOR:  where  P  specifics  a  page  in  the 

flow  diagram,  and  L  designates 
a  statement  number  in  the  program 
listing  or  a  reference  point  in  the 
flow  diagram 

CONNECTOR:  where  X  implies  a  continuation 
of  the  diagram  to  the  next  page 


INPUT/OUTPUT  OPERATION 


PUNCHED  t 


DISK 
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APPENDIX  A 


SUBROUTINE  ALREAD  PROGRAM  LISTING 

SUBFOUTTNE  ALR  E  A  D  (T  ST  A  RT,  TSTOP  ,T  LIFT ,  I  NT  ARC ,  I NPAT ,  IPOL  ,  NO  P  H  A,  NPTS , 
1  DFPG , NEW  P  AS , NRG,  ISTGAT) 

Did  EN  SION  ALT  (5)  ,ABT  (3,5)  ,  A  2  (6)  ,CALCON(  19)  ,  DPPG  (2)  ,  EL  ( 6)  ,  IHD  (1  3 )  , 
IlflT  (5)  ,  INGATE  (5, 3)  , IPRIOR(5)  ,ISENS  (6)  ,ISLIDE(5,3)  ,  ISP  AC  (5,  3)  , 
2TTARDT  (5, 3, 3)  ,  IT  ARC  (5)  ,  TWAVE  (6)  ,LOC(5,  3,4)  ,  NODE  (5,  3)  ,  HEX  (2)  , 

3 N BODES  (5,3)  ,NPTEST  (2)  ,NSAMP(5,3)  ,NSABPT(5)  ,POL(4)  ,  POWER  (6)  , 

4RAD (6,5)  , R  ANGE (5)  ,TYPE(2)  ,VEL(5)  ,IPPAR(22) 

COK1CN/TR  E  AE/LN ,  IELG  ,  IBTRHC  ,  FUTR  HD  , FBTRBA, FBT  RT  G,  FMTR  BI,  F1TRSP, 
1FNXSEC,FB ARSG,FBRR 11, FBGLOT, EBCHAF, FBBSBC, PBASLP, FBABP  (6)  ,FBPHA (6) 
2,NANE(25)  ,NT  (24)  ,IX{24)  ,TABP(128,6)  , TPH  ( 1  28,  6)  ,  ITEN  (2000) 

conncN/TTncon/iHns,iFs 

COnnCN/RDCOBT/TIBES  (300 ) , XSPHA  ( 1 20,  30  0) , RANGKB (300)  ,ALSAV  (30  0)  , 

1 IRGA  (120)  , N P PG 

EQUIVALENCE  (IHD  (1)  ,IDREC)  ,  (IH  D  (2)  ,  LREC)  ,  (IHD  (3)  ,  ITGTB1)  , 

1  (IHD  (4)  , TTGT  B  2)  ,  ( I H C  (5 )  ,  IF  PC) ,  (I  HD (6) , I  PR  I )  ,  (IH  D (7)  ,  N BLR D)  , 

2  (IHD  (9)  , BACYBA)  ,  (I HE  (9)  ,NTARG)  ,  (IHD  (10)  ,  NTDBA )  ,  (IHD  ( 1  1)  ,  RBI  NOR)  , 

3  (IHD  (12)  ,  BICYBA)  ,  (IHD  (13)  ,LBICY)  ,  ( I FP AR  ( 1)  ,  FBT R HD)  , 

4  (PAD  (1,  1)  ,  PO  WF  R  { 1 )  ),(RAD(1,2),AZ(1)),  (RAD  (1 , 3)  ,  EL  ( 1)  )  , 

5  (TTARDT  (1  ,  1,  1)  ,BODE  (1,  1)  )  ,  (ITARDT(  1,  1, 2)  ,ISPAC  (1  ,  1)  )  , 

6  (TTARDT  (1  ,1,3)  ,  NS  A  BP  (1  ,  1)  )  ,  {T  JIT  (  1)  ,  IGCHG)  ,  ( IB T  ( 2)  ,  IGA  IN)  , 

7  (IBT  (  3)  ,  NPLS)  ,  (IBT  (4)  ,  I  SEB  A )  ,  (I  HT  (5  )  ,  NBP)  ,  (ABT  (1,  1)  ,  RANGE  (  1)  )  , 

8  (ABT  (1 ,2)  ,VEL(1) )  ,  (ABT  (1, 3)  , ALT  (1) ) 

DATA  PCON,RKB,VKB/1.0E6, 1.87370 3E- 3,4.4672E-4/ 

DATA  POI/'LC  RC  AZ  EL  •/ 

LATA  TYPE/' VHF-UHE-'/ 

DATA  NET, NPTEST/0, 1,300, 150/ 

DOUBLF  PRECISION  DR  A NG , FI  NT  I B, F PG , G BT I  BE , PRF,T I B E (6)  , T IB ES , T IBOL D, 
1TLTFT,TSTART,TSTOP 
INTEGFR*  2  ITEfl 
TF(NFWOAS.GT. 2)GO  TO  2000 
IF  (  N PG . GT . 120) NRG  = 1 20 
0"TIBE=TSTART-1.0 
N  PT  S  =  0 
TTNC  =  0 . 0 
NT  E  R  R  =  0 

CALL  THEAC  (NEWPAS, 69960) 

DO  60  1-1,6 

IF(IFPAR  (I)  .GT.O) GO  TO  60 
WRITE  (6,4  0)  NABE(I) 

40  PORfl  AT ( ’  FORBAT  TABLE  ',A4,'  WAS  NOT  FOO ND  -  RUN  ABORTED.') 

NFW  P AS=  5  5 
RETUEN 

60  CONTINUE 

80  DO  ICO  1=2,19 

CAT  CON  (I) =GET  (FBXSFC, IBTRHD, I) 

100  CONTINUE 
TLIFT=0. 0 

ITL= IGET  (FBGLOT, IBTRHD,  1) 

TF(ITL.NE.2) GO  TO  140 
I HBS  =  IG FT  (FBGLOT, IBTRHD, 2) 

TES=  IGET  (FBGLOT,  IBTRHD,  3) 

CALL  GBTUPK  (TUFT) 

GO  TO  140 

120  CALI,  BREAD  ( 1) 

140  IHD  (1)  =  IGET  (FBTRHD,IDTRHD,  1) 
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IF  (  (IDREC.LT.  127)  .AND.  (IFLG.  HE.  3)  )  GO  TO  180 
WRITE  (6 ,  160) GMTIME 

160  FORMAT ( •  END  OF  FItE  H  ^COUNTERED  OK  INPUT  TAPE  AFTER  TIME  ',?12„4) 
NFWPAS=44 
RETURN 

180  IF  (ICREC.  NE.  2)  0,0  TO  120 
TIM0LD=GMTIl1F 
DO  2C0  1=2,13 

THD  (I)  =  IGET{FMTRHD,  IBTRHD,  I) 

200  CONTINUE 

IHMS=IGET(FMTRHD,IBTRHD,14) 

IFS=IGET (FMTRHD, IBTRHD, 16) 

CALL  GMTUPK  (GMTIME) 

IF(GMTIME.T T.  (TSTART-0.20) ) GC  TO  120 

IF  (  (GMTIME.  GT.TIMOLD)  .AND.  (M  AC  Y  B  A .  EQ .  0)  .AND.  (NTDBA.EQ.  0)  )  GO  TO  120 
NCONT=0 

IF  (  (  (DADS  (TIMOLD-G  MTI  ME)  )  .GT.0.01)  .AND.  (NT  DBA.  NE.  0)>  GC  TO  240 
IF ( NTERR . NE. 0) GO  TO  120 
NCO  NT= 1 
GO  TO  1200 

240  NCHAN=IFPG*IPOL 

IF  (I FPG . EQ. 1) NCH  AN  =  NCH AN ♦ 3 

NFPG=TFPG* 1 

DFPG  (1)  =TYPE  (N  FPG) 

DO  280  1=1, NOPHA 
INDEX=  (I-  1)  *6*  NCH  AN*  12 
IF(IFPAR  (INDEX)  .NE.O) GO  TO  280 
WRITE  (6,40) NAHE(INDEX) 

N  EWP  A  S  =  66 
RETURN 

280  CONTINUE 

FINTIM=GMTIME 
TIMOI.D=GMTIME 
DO  300  1=1,5 
NSAMPT  (I) =0 
300  CONTINUE 

PPP=FCON/FLOAT (IPRI) 

TINK=1 ,/PRF 

IF  (NBLRD.EQ.O) GO  TO  600 
DO  600  1=  1 , NBL RD 
IBTRMA=IBTRHD*MACYBA*24*(I-1) 

DO  4C0  K=1, 3 

RAD  (I,K)=GET(FMTRMA,IBTRMA,K) 

400  CONTINUE 

I8MS  =  IGET  (FMTRBA,IBTRMA,4) 

I FS= IG ET  (FMTRMA,  IBTRMA,6) 

CALL  CMTUPK  (TTME (I)  ) 

ISENS  (I) =IGFT(FBTRMA,IBTRMA,6) 

IVAVE(I) *IGET (FMTRMA, IETRMA, 7) 

IF(TWAVE(I)  .  EQ.4)  I  WAVE  (I)  =3 
IF  (IWAVE  (I)  .GT.  3)  TWAVE  (I)  =0 
600  CONTINUE 

KONGET=  3*  (NCHAN-1)  ♦IWAVE  (1) *1 

CALADD=CALCON  (KONGET)  -POWER(  1)  -  FLOAT  (  10*ISENS(1)  ) 

600  IF  (NTARG.EQ.O)GO  TO  120 
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ITP IK=0 
NBPTST=0 

DO  900  T= 1 , NTA  RG 
IBTRTG=IBTRHD*NTDBA*12*  (1-1) 

TTARG  (I) =IGET (FflTRTG,  I BTRTG, 1) 

IPRIOR (I) =IGET (FflT  RTG ,  I  ETRT  G,  2) 

N  K  =  1 

DO  800  K  =  1 , 3 
I  K=  N  It  ♦  1 
N  K  =  IK+  2 
DO  7C0  L=  1 , 3 

ITARDT  (I ,  K, L) = IGET ( FflT  RTG, I BTRTG, L*  IK) 

700  CONTINUE 

IF  ( I SPAC  (I ,  K)  . EQ.254)  ISPAC  (I  ,K)  =-1 
ISP  AC (I , K) -2** (ISPAC (I, K ) ♦HIX (N  PPG) ) 

TF  (ISPAC  (I,K)  .EQ.0)ISPAC(I,K)  =  1 
N  SA  APT  (I)  =NSAflPT  (I)  *  NS  A  HP  (I  ,K) 

NflODES  (I,  K)  =0 
DO  780  1=1,4 
LOC  (I,K ,L) =0 
flODUfl=flOD  (NODE  (I, K)  ,2) 

IF (MO DO fl . F 3 . 0 ) GO  TO  760 
NflODES  (I  ,  K)  =  NflODE  S  ( I ,  K)  *1 
LCC (I,K,L)=NflODES ( I  ,K) 

760  NODE  (I, K)=HODE  (I, K)/2 
780  CONTINUE 

ISLICE (I,K) =2*NSAHP (I , K)  *NflODES  (I,R) 

NBPTST=NBPTST*ISLIDE (I,K) 

800  CONTINUE 

DO  820  J=  1, 3 

ING  ATE  (T,  J) =IGET (FflTRTG, I BTRTG, J*  1  1) 

820  CONTINUE 

IF(ITARG(I) .EO.INTARG) ITFIK=I 
900  CONTINUE 

N  PPTST=  N  B  PTST- HOD (NBPTST,6) 

N  TERR  =  0 

IF (ITPIK. GT.O) GO  TO  960 

WRITE  (6, 9 20) G ATI  HE, I  NT ARG, (ITARG (I)  ,I=1,NTARG) 

920  FOR fl  AT  ( '  AT  TIRE  =  ',F12.4,»  TARGET  ',12,'  IS  NOT  ON  THE  TAPE,  TA 
1RGETS  AVAILABLE  ARE  ',513) 

IBTRflI=TBTRHD*1TCTBA 
NPP  =  IGET (PflTRflI,IBTRflI ,5) 

IF (NEP. EQ.NEP1ST) GO  TO  10000 
NTEF  F=  1 
GO  TO  120 

960  IF (LCC (ITPIK, INPAT, IPOL)  . GT. 0 j  GO  TO  1020 

WRITE  (6, 1000) 

1000  PORfl  AT ( '  EITHER  POIARI ZA^tcn  OR  PATTERN  CHOSEN  IS  NOT  AVAILABLE') 
GO  TO  10000 

1020  IF  fl ST GAT. GT.NSAflP( ITPIK, IN  PAT)  )  ISTGAT=1 
IF  (ISTGAT.LT. 1)ISTGAT=1 
N  STEST=NSAflP ( ITPI K , I NP AT) -I STG AT* 1 
NSTO AT= ISTGAT 
TPAT=TNP AT 
1  =  1 
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10U0  IRGA  (I)  *  INGATE ( ITPIK,  IPAT)  ♦ ( NSTG AT- 1 ) * ISPAC (ITPI K , IPA T) 

NSTG  AT=  1 

1060  IF(I.GE  .NRG) GO  TO  1100 

IF (  *1+1)  .GT.NSTEST)  GO  TO  1080 
1=1  +  1 

IRGA  (I )  =1 RG A  (T-1)  +ISPAC  (ITPIK, IPAT) 

GO  TO  1C60 
1080  I PAT  =  TP  AT ♦ 1 

IF (IPAT.GT. 3)GO  TO  1100 
IF(NSANP(ITPIK,IPAT) . LE. 0) GO  TO  1080 
NSTEST=NSTEST+NSANP (ITPIK,  IPAT) 

T  =  I  ♦  1 

IF  (I.  IE.  NRG)  GO  TO  1040 
1  =  1-  1 

1100  I  PA  T  =  T N  P A T 
N  RC>=  I 

DFPG  (2)  =POL  (TPOi) 

1200  IBTRHI=IBTRHD+HICYEA 
DO  2200  HTN= 1 , NHINOF 
IF(HIN.LE.  1)  GO  TO  1220 
IBTRHI=IBTRNI*LHICY+NPLS*NBP 
1220  DO  1 2*40  1  =  1,5 

INT  (I)  =  IGET  (FHTRHI  ,IBTRHI  ,1) 

1240  CONTINUE 

IF (NBP.  EQ. NBPTST) GO  TO  1290 
WRITE  (6,  1260)GNTIHE,HIN,NBP,NBPTST 

1260  PORN  AT ( '  AT  TINE  =  ',F10.4,',  HINOR  CYCLE  ',11,'  SOHETHING  IS  SCRE 
1WFD  (IF  NBP  =  ',15,'  IT  SHOULD  BE  ',14,'  -  SHOOT  H A RTOGEN SIS  • ) 

GO  TO  1330 

1280  IF (NPLS.GR.  1) GO  TO  1340 

WRITE  (6, 1  320)  NPLS, NIN, NNINO  B 

1  320  POR  NAT  ( '  NPLS  =  • , I  2, 3X, • NINOR  CYCLE  ', 11,31,11,'  HINOR  CYCLES') 
1330  GNTINE=GNTINE+TINC 
GO  TO  120 

1  3 <4 0  TINC  =  FLOAT  (NPLSJ/PRF 
DO  1500  L=1 , NTARG 
IBTRNT=IETRNI+ ( L— 1 ) *12 
DO  1<400  T  =  1 , 3 

ANT  (L,I)  =GET(FNTRHI,IBTRNT,I+5) 

1400  CONTINUE 

VFL  (I) = VEL  (L) *VKN 
1800  CONTINUE 

IBTRHT=IBTRNI  +  (ITPIK- 1)  *12 
ISDREl=IGET  (FNTRNI,  IBTRHT,  R) 

NSDUH=NSANPT(ITPIK) 

IF(NRG.LT.NSCUN) NSDUH=NRG 
DRANG= RANGE (ITPIK) 

IF (DRANG. GT.O. 0) GO  TO  I860 
WRITE (6, 1540) GNTINE,D RANG 

1540  FORNAT  ( '  AT  TINE  =  ',F15.4,»  THE  RANGE  =  ',1PE20.6) 

DPANG=1 . 0 

1560  DPANG=4C. *DLOG 10 (DRANG) 

CALI B= DRANG* FLO AT (IGA  IN)  +CALADD 
DO  2000  1=1, NPLS 

IF( (NIN.GT.1) .OR.  (I.GT.1))  FI NTI H=FI NTIN  +  TINK 
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I F (TST  A  RT . GT. PINT  I M) GO  TO  2000 
NFTS=NPTS  + 1 

IF (NPTS.GT. 1) GO  TO  1600 
1580  TIMES(NPTS) =FINTIM 
GO  TO  1620 

1600  IF  (  (I.  EQ.  1)  .  AND.  (MTN.EQ.  1)  .  AND.  (NCONT.  EQ.  0) )  GO  TO  1580 
TIMES (NPTS) =TI MES  (NPTS-1) *TI UK 

16  20  RANG KM  (NPTS)  =  R  ANG  J  (ITPIK)  *RKM*VEL( ITPIK)  *  (TIMES  (NPTS)  -FINTIM) 
AI.S  A  V  (NPTS)  =»LT  (ITPIK) 

ISAMPT=IBTRHD  +  ISDRM  ISDREU  (1-1)  ♦N BP 
NSTEST  =  N SAMP  (ITPIK , INP AT)  - 1 STGA T* 1 
I  AD=0 

IF(IN£AT.EQ. 1)GO  TO  1680 
JST=?NP AT-1 
DO  1660  J=1,JST 
IAD=IAD+I SLIDE  (ITPIK,  J) 

1660  CONTINUE 
1680  L=I STGAT -1 

DO  1800  K -  1 , NRG 
L=L  +  1 

IF (K. LE. NSTEST) GO  TO  1700 
IAD-IAD+ISLIDE (ITPIK, IPAT) 

IPAT=1PAT*1 

NSTEST= NSTEST »NS AMP (ITPIK, IPAT) 

L=  1 

1700  IPIK=ISAMPT+I AD+2*  ( IOC  ( IT  PI  K  ,1 P  AT ,  I  POL)  -  UNMODES  (ITPIK,  IPAT) 

i*a-i)) 

IAMF=IGET  (FMTFS  P,  IP  IK ,  1) 

IF( (IAMP. LT.  1) .OH.  (IAMP.GT.  128) )  IAMP-1 
1720  XSPHA (K,NP1S) =TAMP (IAMP, KCH AN) ♦ CALI B 
1740  IF (NOPHA. HO. 1) GO  TO  1800 
IPHA-IGET  (F MTRSP, I PIK , 2) 

IF(  (TPHA.GE.O)  .AND.  (IPHA.LE.  127)  )GO  TO  1760 
XSPHA (K, NPTS+150) =0.0 
GO  TO  1800 

1760  XSPHA  (K,NPTS+  150)  =TPH  (IPHAH  ,  NCH  AN) 

1800  CONTINUE 

IPAT=INPAT 

IF(TSTOP.LT.TIMES(NPTS)  )GO  TO  10000 
IF (NFTS. LT. NPTEST (NOPHA) ) GO  TO  2000 
NEWP AS=  99 
RETURN 

2000  CONTINUE 
2200  CONTINUE 
GO  TO  120 

8960  WRITE (6,9980) 

9980  FORMAT  ( •  THEAD  HAS  DEFAULTED  -  RUN  ABORTED.') 

1000J  NEWPAS=0 
RETURN 
END 
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APPENDIX  B 

SUBROUTINE  aLREAD  FLOW  DIAGRAM 
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APPENDIX  C 


-SUBROUTINE  THEAD  PROGRAM  LISTING 

SUBROUTINE  THEAD  (NEWPAS,  *)  VERSION  3/15/71 

INTEGER*2  TTEH 

COHHON/TRE  AC/LN,IFLG,  IADD, ICATFB  (24)  f  NABE  (25)  ,NI(24)  ,  I X  ( 2  4)  , 

IT ABP  (128,6)  ,TPH  (128,6)  ,ITEB(  20  0  0) 

DIHENSION  IBUF  (2048,2)  ,IBUF1  (2048)  ,IBUF2  (2048)  ,IBARR(8) , N  AH  EX (25) , 
1TEH  ( 1) 

EQUIVALENCE  (IBUF1  (1 )  ,  I  BOP  (1 ,  t )  )  ,  ( IBUF2  ( 1)  ,  IBUF  ( 1,  2)  ) 
EQUIVALENCE(TTEH(1)  ,TEH  (1)  ) 

DATA  NAHEX/'TRHD*  ,  '  T  R  B  A  ’  ,  •TRTG*  ,  'TRBI  '  ,  «TRSP'  ,  •  XSEC  ,  •  ARSG'  ,  •  RR1  1  • 
1, 'GLOT* , ' CH  AF ' , ' BSBC*  ,' ASLP'  ,' A  HP1 '  ,  '  ABP2  ' , ' ABP3 ' ,  ' ABP4 ' ,  '  ABP5 ' , 

2 ' AH  P6 ' , • PHA1 ' , • PH A  2* , ' PHA3' , ' PHA4' , •PHA5’ , • PHA6' , 'HDBD'/ 

DATA  BAX/8192/, IT/1/ 

DATA  ICBS/'CBS  •/ 

IF  (NEWPAS. NE.O) GO  TO  40 
CALL  WHICHV  (ID) 

IOUT=6 

TF  (IE. EO. ICBS)  IOUT=8 

NREC-=0 

LN  =  1 

DO  20  1=1,24 
ICATFB  (I)=0 
NABE  (I) =NABEX (I) 

NT  (I)=0 
IX  (I)=0 

20  CONTINUE 

ICATFB ( 10 ) = 1 

NAB  E (25)  =  NABEX  (25) 

40  CALL  BREADS (LN, IBUF 1, I BDF2, BAX, IFL, INDX, LEN,  IFLG , IADD) 

60  CALL  BREAD(LN) 

IF(NEWPAS.NE.O)  RETURN 
NREC=NREC+1 

IF (IFLG .EQ.2) GO  TO  380 
IF(IFLG.EQ.3.AND.IT.EQ.1)G0  TO  60 
IF(IFLG.EC.3)GO  TO  480 

IF((IBUF(1, INDX) /16777216) .EQ.2) GO  TO  4B0 
70  CALL  NAHET(IADD,NABID) 

IF(NABED. EQ.NABEX (25))  GO  TO  100 
WPITE  (10 UT , 80) NABED 
80  FORB  AT ( 1 X  ,A4 , ’  FOUND') 

100  DO  120  1=1,25 

IF (NABED. EQ.NABEX(I)  )  GO  TO  140 
120  CONTINUE 

GO  TO  60 

140  IF (NABED. EQ.NABEX (10)  ) GO  TO  60 

CALL  FORB  (IADD,ITEB (IT)  ,1 B , NABED, NT EB , 6260) 

IF(NABED. EQ.NABEX (25) ) GO  TO  420 
DO  160  1=1,12 

TF(NABE(I) .EQ. NABED) GO  TO  220 
160  CONTINUE 

DO  180  J=  13 , 1 8 
I=J- 12 

IF (NABE (J) . EQ. NABFD) GO  TO  300 
180  CONTINUE 

DO  200  J=19,24 
T= J-  18 


IF  (NAME  (J)  .F,Q.  NAMED)  GO  TC  340 
200  CONTINUE 
GO  TO  60 

220  ICATFM  (I) =IB 

WRITE (I0UT,240) NAMED , NTEM 

240  FORMAT  {'  FORM AT= ' , A4, '  STORFD  IN  COMMON  NTEN=',I4) 

IT  (I) =  IT 
NI  (I)  =NTEM 
IT= IT+6  *NTEM 
GO  TO  60 

260  WRITF(6f 280) NAMED, NTEM 

280  FORMAT  ( •  NAME  =  'A4,'  NTEM  ='I5,'  *ERFOS*  FORMAT  TABLE  LIMITED  T 
10  700  ITEMS  CR  FORMAT  TABLE  HAS  0  LENGTH') 

GO  TO  4  80 
100  JT=IT 

ICATFM  (J) =IB 
DO  320  J=  1 ,  127 
JT= JT+6 

TAMP  (J,I)=TEH(  (JT*5)/2) 

320  CONTINUE 
GO  TO  60 
340  JT=  IT 

ICATFM  (J)=IB 
DO  160  J=  1,  128 
JT=  JT+6 

TPH  (J,I )  =TEN  ( (JT*6)  /2) 

160  CONTINUE 
GO  TC  6 C 

180  WRITE  (IOUT, 400) NFEC 

400  FORMAT  (•  PARITY  ERROR  READING  FORMAT  RECORD', 16) 

GO  TO  60 

420  C’iLL  HDRR(IADD,ITTP,IDAFR) 

WMTE  (IOUT,  440)  ITYP 
440  FORMAT ( '  TYPE  • ,12) 

IF ( (ITYP. EQ. 1)  .OR.  (ITYP. EQ. 2) ) CO  TO  60 
WRITE  (6,460) 

460  FORMAT  ( '  TAPE  NOT  TRANSCRIPTION  TAPE  JOB  TERMINATED  BY  THBAD' ) 
RETURN  1 

480  DO  6 20  1=1,24 

IF ( ICATFM  (I) .EQ.O) WRITE (I0UT,500) NAME! (I) 

SCO  FORM  AT( IX, A4,'  NOT  FOJND') 

620  CONTINUE 

WRITE (6,540) 

540  FORMAT('  THRAD  CONPIETE') 

RETURN 

END 
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APPENDIX  E 

SUBROUTINE  THEAD  OUTPUT 
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APPENDIX  F 


SUBROUTINE  GMTUPK  PROGRAM  LISTING 


100 

200 


SUBROUTINE  GMTUPK (TIME) 

DIMENSION  SFCMUL  (6,  2)  ,T  (2)  ,LIM  (2)  ,  ISHIFT(6,2) 
COMMCN/TI NCOM/ITM  (2) 

DOUBLE  PRECISION  TI  ME  ,  T  ,T  SUM 

CA* A  SECMU L/36000. 0, 36 00. 0,6 CO. 0,60. 0,10. 0,1. 0,0. 

10.0,0.0/ 

DATA  I SHIFT/ 1048 576, 65 5 3 6, 4096, 2 56, 16, 1,256,  16,  1, 
DATA  LIM/6,3/ 

DO  200  J=  1 , 2 
TSUM  =  C .  0 


1,0. Cl, 0.001, 0.0, 
1,1,  1/ 


IDIG=  IT  M  (»T) 

ILIM=LI H (0) 

DC  ICO  I  - 1 , ILI M 
NSLIEE=IDIG/ISHIFT  (I,J) 

TSUM=T  SUM  ♦FLOAT  (NSLIDE)  *SECMUL  (I  ,.1) 

IDIG  =  IDIG-NSLIDE*ISHIPT  (I  ,J) 

CONTINUE 

T  (J)  =TS  0  M 

CONTINUE 

TIN  E=T  (  1)  +T  (2) 

RETURN 

END 


